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Abstract 

We present a brief review on the Raychaudhuri equations. Beginning with a summary of the 
essential features of the original article by Raychaudhuri and subsequent work of numerous authors, 
we move on to a discussion of the equations in the context of alternate non-Riemannian spacetimes 
as well as other theories of gravity, with a special mention on the equations in spacetimes with 
torsion (Einstein-Cartan-Sciama-Kibble theory). Finally, we give an overview of some recent 
applications of these equations in General Relativity, Quantum Field Theory, String Theory and the 
theory of relativisitic membranes. We conclude with a summary and provide our own perspectives 
on directions of future research. 



I. THE BEGINNINGS 



About half a century ago, General Relativity (GR) was young (just forty years old!), 
and even the understanding of the simplest solution, the Schwarzschild, was incomplete. 
Cosmology was virtually in its infancy, despite the fact that the Friedmann-Lemaitre- 
Robert son- Walker (FLRW) solutions had been around for quite a while. The question 
about the then-known exact solutions of GR, which worried the serious relativist quite a 
bit, concerned their singular nature. Both the Schwarzschild and the cosmological solutions 
were singular. It is well-known that the creator of GR, Einstein himself, was quite worried 
about the appearance of singularities in his theory. Was there a way out? Was it correct to 
believe in a theory which had singular solutions? Were singularities inevitable in GR? 

It was during these days in the early 1950 's, Raychaudhuri began examining some of 
these questions in GR. One of his early works during this era involved the construction 
of a non-static solution of the Einstein equations for a cluster of radially moving particles 
in an otherwise empty space A year before, he had also written an article related to 
condensations in an expanding universe where he dealt with cosmological perturbations 
(in a sense, this article deals with what is today known as structure formation) . Subsequent 
to these papers, in 1955, appeared Relativistic Cosmology I [Sj], which contains the derivation 
of the now-famous Raychaudhuri equation. 

Fifty years hence, the Raychaudhuri equations have been discussed and analysed in a 
variety of contexts. Their rise to prominence was largely due to their use (through the notion 
of geodesic focusing) in the proofs of the seminal Hawking-Penrose singularity theorems of 
GR. Today, the importance of this set of equations, as well as their applicabihty in diverse 
scenarios, is a well-known fact. 

This article is a brief review on these equations. We shall deal with some selected aspects 
in greater detail. We, of course, would like to emphasize that there are many topics which 
we leave untouched or, barely touched. We hope to do justice to these in a later, and more 
extensive article. 

The overall plan of this article is as follows. In the remainder of this section, we shall recall 
the basic ideas and results in the 1955 paper and the singularity theorems. The next section 
introduces (with illustrative examples) the kinematical quantities (expansion, rotation and 
shear) which govern the characteristics of geodesic flows and also outlines the derivation and 
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consequences of the equations. Section 3 considers the equations in alternative Riemannian 
and non-Riemannian theories of gravity {R + jSR"^ theory and the Einstein-Cartan-Sciama- 
Kibble (ECSK) theory, in particular). In Section 4, we give a glimpse of the diverse uses of 
these equations in contexts within, as well as outside the realm of GR. Finally, we present 
a summary and provide our perspectives on possible future work. 



A. The original 1955 paper 

The derivation of the Raychaudhuri equation, as presented in the 1955 article, is some- 
what different from the way it is arrived at in standard textbooks today. It must however be 



mentioned, that in a subsequent paper in 1957 {4], Raychaudhuri presented further results 
which bear a similarity with the modern approach to the derivation. Let us now briefly 
summarise the main points of the original derivation of Raychaudhuri. 

(a) Raychaudhuri 's motivation behind this article is almost entirely restricted to cosmology. 
He assumes the fact that the universe is represented by a time-dependent geometry but does 
not assume homogeneity or isotropy at the outset. In fact, one of his aims is to see whether 
non-zero rotation (spin), anisotropy (shear) and/or a cosmological constant can succeed in 
avoiding the intial singularity. 

(b) The entire analysis is carried out in the comoving frame (in the context of cosmological 
line elements) -the frame in which the observer is at rest in the fluid. 

(c) The quantity Rf (spacetime coordinates in the 1955 paper are labeled 

with the fourth one being time), is evaluated in two ways-once using the Einstein equations 
(with a cosmological constant A) and, again, using the geometric definition of R\ in terms 
of the metric and its derivatives. In the second way of writing this quantity, Raychaudhuri 
introduces the definitions of shear and rotation. 

(d) Finally, equating the two ways of writing R\ the equation for the evolution of the 
expansion rate is obtained. Note that the definition of expansion given in this paper refers 
to the special case of a cosmological metric. 

(e) Apart from obtaining the equation, the article also arrives at the focusing theorem 
(though it is not mentioned with this name) and some additional results (the last section). 

In the same year (1955), Heckmann and Schucking j^, while dealing with Newtonian 
cosmology arrived at a set of equations, one of which is the Raychaudhuri equation (in the 
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Newtonian case). Prompted by this work, Raychaudhuri re-derived his equations in a some- 
what different way in an article where he also showed that Heckmann and Schucking's work 
for the Newtonian case could be generalised^without any problems to the fully relativisitic 
scenario. It must also be noted that Komar [6], a year after Raychaudhuri's article appeared, 
obtained conclusions similar to what is presented in Raychaudhuri's article. Raychaudhuri 
pointed this out in a letter published in 1957 [7|. 

Subsequently, in 1961, Jordan, Ehlers, Kundt and Sachs wrote an extensive article on the 
relativistic mechanics of continuous media where the derivation of the evolution equations of 
shear and rotation seem to appear for the first time jsl . Furthermore, for null geodesic flows, 
the kinematical quantities: expansion, rotation and shear (related to the so-called optical 
scalars) and the corresponding Raychaudhuri equations, were first introduced by Sachs [9|. 

The Raychaudhuri equation is sometimes referred to as the Landau-Ray chaudhuri equa- 
tion. It may be worthwhile to point out precisely, the work of Landau, in relation to this 
equation. Landau's contribution appears in his treatise The Classical Theory of Fields [l^ 
and is also discussed in detail in p, 11|. Working in the synchronous (comoving) reference 

dt ' 

y, where 7 is the 

expression for Rq and then, an inequality + | (Xa)^ ^ 0. While deriving the inequality. 
Landau implicitly assumes the Strong Energy Condition (though it is not mentioned with 
this name). Then, using it, he is able to show that 7 must necessarily go to zero within a 
finite time. However, he mentions quite clearly that this does not imply the existence of a 
physical singularity in the sense of curvature. Though Landau's work captures the essence of 
focusing, he does not explicitly mention geodesic focusing-moreover, he does not introduce 
shear and rotation or write down the complete equation for the expansion. 



frame. Landau defines a quantity Xafs = where 'ja/s is the 3-metric. Subsequently, 

using the fact that Xa — where 7 is the determinant of the 3-metric, he writes down an 



Even though it was mentioned in llj and 



Raychaudhuri's contribution found its true 



recognition only after the seminal work of Hawking and Penrose which appeared a decade 
later. It was at that time, along with the proofs of the singularity theorems, the term 
Raychaudhuri equation came into existence in the physics literature. 
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B. The singularity theorems 



It must be mentioned that Raychaudhuri did point out the connection of his equations 
to the existence of singularities in his 1955 article. However, more general results (based on 
global techniques in Lorentzian spacetimes) appeared in the form of singularity theorems fol- 
lowing Penrose's work 12 1 and, then, Hawking's contributions Q, [l^. The crucial element 
of the singularity theorems is that the existence of singularities is proved by using a minimal 
set of assumptions (loosely speaking, these are : Lorentz signature metrics and causality, the 
generic condition on the Riemann tensor components, the existence of trapped surfaces and 
energy conditions on matter). In fact, a precise definition of 'what is a singularity?' first 
appeared in the works of Hawking and Penrose. The notion of geodesic incompleteness and 
its relation to singularities (not necessarily curvature singularities) was also born in their 
work. One should also realise that the focusing of geodesies arrived at by Raychaudhuri 
and discussed in much detail in later articles by other authors could be completely benign 
(irrespective of any actual singularity being present in the manifold). Thus, a singularity 
would always imply focusing of geodesies but focusing alone cannot imply a singularity (also 
pointed out by Landau lOj ) • We refrain from discussing the singularity theorems any further 
here-excellent discussion on glob al aspects in gravitation as well as the Hawking-Penrose 



theorems are available in 



15 



16 



17|. 



II. THE GEOMETRY AND PHYSICS OF THE EQUATIONS 

Let us now review the basic ingredients and the derivation of the equations. First, of 
course, we need to know- what do these equations deal with? In a sentence, one may say 
that they are concerned about the kinematics of flows. Flows are generated by a vector 
field-they are the integral curves of the given vector field. These curves may be geodesic 
or non-geodesic, though the former is more useful in the context of gravity. Thus, a flow 
is a congruence of such curves -each curve may be timelike or null or, in the Euclidean 
case, have tangent vectors with a positive definite norm. One does not, in the context of 
these equations, ask, how the flow is generated. In other words, we are more interested, 
in deriving the kinematic characteristics of such flows. The evolution equations (along the 
flow) of the quantitites that characterise the flow in a given background spacetime, are the 
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Raychaudhuri equations. Historically speaking, it is the equation for one of the quantitites 
(the expansion), which is termed as the Raychaudhuri equation. However, in this article, 
we will refer to the full set of equations as Raychaudhuri equations. 



A. Expansion, rotation, shear 

What quantities characterise a flow? If A denotes the parameter labeling points on the 
curves in the flow, then, in order to characterise the flow, we must have different functions 
of A. In other words, the gradient of the velocity field being a second rank tensor is split 
into three parts : the symmetric traceless part, the antisymmetric part and the trace. These 
define for us the shear, rotation and the expansion of the flow. Specifically, 

"^bVa = (Tab + UJab H ^-/^afe© (l) 

where the symmetric, traceless part, the shear, is defined as aab = | (V^fa + Vafb) — :^^/?-a60, 
the trace, expansion, is = Vaf" and the antisymmetric rotation is given as, Uab = 
\ (Vfefa — Vaffe). n is the dimension of spacetime, and hab = Qab ± VaVb is the projection 
tensor (the plus sign is for timelike curves whereas the minus one is for spacelike ones). Also, 
correspondingly, VaV"" = =Fl- We shall discuss the case of null geodesic congruences briefly, 
later. 

The geometric meaning of these quantities is shown through Figure 1 and Figure 2. 
The expansion, rotation and shear are related to the geometry of the cross sectional area 
(enclosing a fixed number of geodesies) orthogonal to the flow lines (Figure 1). As one moves 
from one point to another, along the flow, the shape of this area changes. It still includes 
the same set of geodesies in the bundle but may be isotropically smaller (or larger), sheared 
or twisted. The analogy with elastic deformations or fluid flow is, usually, a good visual 
aid for understanding the change in the geometry of this area. A recent, nice discussion is 
available in Earlier references where these quantities are explained in quite some detail 
are 19|,|2oi- 
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Area enclosing a set of flow lines 




FIG. 1: The cross sectional area enclosing a congruence of geodesies 



Isotropic Expansion 



Area at X j 




Area at X 



Area atX . 



Rotation 




Area at X . 



Shear 




Area at Area at X, 2 



Note similarity with 

(i) Elastic deformations 

(ii) Fluid flow 



FIG. 2: Illustrating expansion, rotation and shear 



B. Examples 



It is useful to illustrate these quantities with a set of examples. We first choose to work 
with Schwarzschild spacetime. Our examples here will involve (i) rotation-free timelike 
geodesic flows and (ii) timelike flows with all three kinematical quantities non-zero. We 
focus on examples with non-zero shear and rotation because these are not usually available 
in standard texts on GR. Our choice of examples are primarily based on the problems 
suggested in a recent monograph by Poisson 3] • In a third example (iii) we discuss briefly 
a timelike geodesic flow in the FLRW universe. Finally, in (iv) we briefly deal with geodesic 
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flows in wormhole spacetimes. 
(i) From the Frobenius theorem [18|, we know that hypersurface orthogonal vector fields 
must necessarily have zero rotation though the shear can have non-zero components. In 
Schwarzschild spacetime, we construct a congruence which has the above properties (i.e. it 
is irrotational but has nonzero shear and expansion). Consider the vector field: 



u^d. = ^^d,±f-fdr (2) 

r 

The geodesies corresponding to the above vector field are marginally bound (i.e. Ut = — 1) 
and the upper and lower signs refer to outgoing and incoming geodesies. It is easy to show 
that the above vector field can be written as Ua = da(f), where 0(x°) = constant would 
represent the hypersurface with respect to which u"' is orthogonal. In fact, (p is the same as 
the new time T (the proper time as measured by a freely falling observer starting from rest 
at infinity and moving radially inward) used in the Painleve-GuUstrand representation [l^ 
of the Schwarzschild line element. 

It is easy to calculate the expansion, which turns out to be: 



, 3 2M 

e = ±2 V — 

Notice that the expansion is positive for outgoing and negative for incoming geodesies. We 
can also find the nonzero shear tensor components which are given as: 



2M 2M 



<ytt = T^r\ ; crrr = T- ,2 

\ r fi _ 2M_Y 



2M 

r-2 



^tr - ^ _ 2M ; ~ ^ V ~ ~ sin^^'^'^ ^ ' 

r 

(5) 

One can then check that the Raychaudhuri equations (given in the next sub-section) hold 
with the above expressions for the shear and expansion. 

(ii) We now move on to an example where all three kinematical quantities are non-zero. 
Consider the following vector field in Schwarzschild spacetime: 



^""da^^—f d, + J-,dA (6) 
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where M is the usual mass. It can be verified that the geodesies corresponding to this 
vector field are timehke and they are circular (r is a constant). We intend to calculate the 
expansion, rotation and shear of the above vector field. The expansion is given as : 



e = cot^ 



M/r3 



3M 
r 



(7) 



Notice that the expansion is positive in the northern hemisphere and negative in the southern 

hemisphere. 

Following the definition, one can show that the rotation tensor for this vector field is 
given as: 



M 1 



6M 



OJtr — 



47^2 



^1 _ 3M 



(8) 



and the shear tensor is: 




One can further verify that here too, the Raychaudhuri equations (given below) are satisfied 
for the above quantities. One can also evaluate — u"^ and show that it is positive for 
r > 3M. The expansion is defined for this domain of r (> 3M) and can diverge to negative 
infinity (focusing) at ^ = tt (south pole). 

(iii) As a third example, we now quickly discuss the expansion, rotation and shear with 
respect to the vector field u'^da — dt in the standard cosmological fine element (FLRW). The 
shear and rotation are identically zero. The expansion is given as: 

e = 3^ = ^^(V^) (10) 

Note that is the volume of the expanding 3-space. Hence the term 'volume' expansion 
is also used in the literature. Raychaudhuri, in his original article, defined the expansion 
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in the this way. However, his treatment did include nonzero shear and rotation because he 
did not assume, to start with, maximally symmetric metrics on spatial slices representing 
R^, or H^. It may be mentioned here that the equation for the expansion reduces to the 
equation for - = (p + 3p). 



(iv) Our final example concerns the case of geodesic flows in a traversable wormhole |21|. 
It is known that traversable wormholes require energy-condition-violating matter. These 
are non-singular spacetimes where the spatial slices resemble two asymptotically fiat regions 
connected by a throat. Thus a geodesic congruence passing through the throat from one 
asymptotic region to the other would necessarily tend to focus first, not quite reach a focal 
point, and then defocus. For a typical wormhole, the line element is given as: 

ds^ = -x^{l)df + dl^ + r2(/) (dO^ + sin^ ed(P^) (11) 

where, for a wormhole, x(0 is nonzero and finite for all I and r{l = 0) = 60 with r{l — »• 
±00) ~ /. It is easy to check that the expansion is proportional to (the prime denoting 
a derivative w.r.t. 1). Thus, the expansion never becomes negative infinity and thus there 



is no focusing. One can work out the expansion for a typical example using r{l) = y&o + 
(the so-called Ellis wormhole) j21|. 

C. The equations and the focusing theorem 

We now turn towards writing down the evolution equations for the expansion, shear and 
rotation along the flow representing a timelike geodesic congruence. A fact worth mentioning 
here is that, these evolution equations (and their generalisations) are essentially geometric 
statements and are independent of any reference to the Einstein field equations. 



The modern (textbook) way to derive these equations (see 15|) is as follows. Consider 
the quantity v'^VcBab (where Bab = V^fa). Evaluate this as an identity and then split it 
into its trace, antisymmetric and symmetric traceless parts. The equations that emerge are 
the ones given below (for n = the dimension of spacetime = 4). 

^ + W + ^'-^' = -RabV^^' (12) 

ctA 3 
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da, 



ab 



dX 



+CcbadV%'^ + ^Rab 



(13) 



du^ 



ab 



-ecu, 



'^Cr''[b^a]c 



ise that these are not equations 
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2J] we find the usage Ray- 



2 'ab - ^f,UJa]c (14) 

where cr^ = CafeO""*, cu^ = UabUj""^, Ccbad is the Weyl tensor and the quantity Rab = hachbdR'^'^ ■ 
\habhcdR'^'^ ■ 

There few points to note here. Firstly, one must real: 
but, essentially, identitites. Hence, in some references 22, 
chaudhuri identity or Codazzi-Raychaudhuri identity (in the context of surface congruences 
to be discussed later) which is, indeed, rigorously correct. The identities, however become 
equations once we use the Einstein equations or any other geometric property (e.g. Einstein 
space, or vacuum, etc.) as an extra input. However, we shall continue to use the term 
equations in this article. 

Furthermore, the equations are coupled, nonlinear and first order. The equation for the 
expansion is of central interest (in the context of the singularity theorems) and it is rather 
straightforward to analyse. In mathematical parlance, it (the equation for the expansion) 
is known as a Riccati equation. Such equations can be transformed into a second order 
linear form (more precisely a Hill-type equation or a harmonic oscillator equation with a 



time- varying frequency) 25|, l26|. Redefining Q = 3yr one gets: 



d^F 



+ \ {RabV^v" + 







(15) 



■ 3 

The analysis of the expansion equation can be done using the above form. One notes that 
the expansion B is nothing but the rate of change of the cross-sectional area orthogonal 
to the bundle of geodesies. Therefore, the expansion approaching negative infinity implies 
a convergence of the bundle, whereas a value of positive infinity would imply a complete 
divergence. What are the conditions for convergence? Firstly, for convergence we must have 
an initially negative expansion. Finally, with F' negative we must end up at a zero of F 
(at a finite A), in order to have a negatively infinite expansion 



25 



26j. Thus, the criterion 



for the existence of zeros in F at finite values of the affine parameter is what is required for 
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convergence. Using the well-known Sturm comparison theorems in the theory of differential 
equations one can show that convergence occurs if : 



RabV^'v^ + a^-uj^>0 (16) 

Thus, rotation defies convergence, while shear assists it. The equation for the evolution 
of the rotation Uab, has a trivial solution Uab = 0. The criterion for convergence then 
becomes particularly simple for such hypersurface orthogonal congruences (zero rotation) : 
Rabv"'v^ > 0. This leads to geodesic focusing. 

If we make use of the Einstein field equations and rewrite the Ricci tensor in terms of the 
energy-momentum tensor Rab = Tab — \gabT the the so-called timelike convergence condition 
becomes a condition on matter stress energy. This, given as, {Tab — \gabT^ 'v'^'v^ > is known 
as the Strong Energy Condition (SEC). For a diagonal Tab (with Tqo = p,Taa = Pa) "we must 
have p + Pa> 0, P + J2aPa > if the SEC is to be obeyed. In other words, geodesic focusing 
encodes the simple statement that if matter is attractive, geodesies must be eventually drawn 



towards each other. This seeming 



In the late seventies, Tipler 
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y trivial statement is proved via the focusing theorem. 



26( 1 realised that the assumption of the SEC imposed to 
prove focusing and hence the existence of singularities could be further weakened. Among 
other results, he was able to show how, in the proof of the Hawking-Penrose theorem one 
might replace SEC by the Weak Energy Condition (WEC: Tabv"'v^ > V non-spacelike w"). 
Tipler also introduced in his article, for the first time, the notion of an averaged energy 
condition (the Averaged Strong and Weak Energy Conditions (ASEC and AWEC) which 
are global in nature. For instance the AWEC is obtained by integrating the WEC along a 
non-spacelike geodesic and gives a number {J^^ TabV'^v^dX > 0). The question of whether a 
violation of the energy conditions could lead to a non-singular solution was also addressed 



by him. One should mention here that a few years before Tipler's work, Bekenstein 



27| and 



Murphy [28| had proposed certain non-singular spacetimes. These models turned out to be 
singular, following the above-mentioned generalisations of the singularity theorems. Tipler 
also quotes, in his paper, the well-known Epstein-Glaser-Yaffe theorem in quantum field 
theory where it is said that there can exist a quantum state w.r.t. which the expectation 



value of the stress-energy tensor can be negative [29|. However, he does mention that 
the existence of one state w.r.t. which the (Too) is less than zero cannot really lead to the 
prevention of singularities. Of late, however, energy-condition violations have been discussed 
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in great detail in the context of wormholes [2l| , dark energy [30| , braneworld models [31| and 
semi-classical gravity (quantum field theory in curved spacetimes) 32| . The experimentally 
observed Casimir effect (the feeble attraction between parallel, conducting capacitor plates) 



331] has been cited as an example of the existence of negative energy density (though, truly 



speaking this effect is concerned with negative pressures as opposed to actual negative energy 
densities) . 

The quest for further weakening the criteria for geodesic focusing and thereby making the 
sing ularity theorems stronger was continued later through the work of Borde 3J] and Roman 



35l | . Further work on issues related to geodesic focusing, singularites, energy conditions and 



causality violations have been carried out in |36|, 



37 



38 



39|,|40|. 



D. Null geodesic congruences 

The Raychaudhuri equations for null geodesic congruences were first derived by Sachs [91] 
in 1961. Let us briefly recall the salient features of these equations. It must be mentioned, 
however, that these equations are not very different (in structure as well as consequences 
thereof) from the equations for timelike congruences. 

The central issue in the case for null geodesic congruences is the construction of the 
transverse parts of the deviation vector and the spacetime metric. Assuming an affine 
parametrisation in the sense dx°' = k°'d\ with k"'ka = and k"-C,a = (with C,"" being the 
deviation vector) we realise that we are in trouble because of the above two orthogonality 
relations. Naively writing hab = gab + kak^ will not work here {k"'hab 7^ 0). The transverse 
metric is thus constructed by introducing an auxiliary null vector A^*^ with kaN"" = — 1. 
(the choice of minus one is by convention, the essence is that the quantity must be non- 
zero). If we choose ka = —daU {u = t — x) then we can have A^^ = —\daV and hence 
hab = dab + kaNb + kbNa- This satisfics k°-hab = and N^hab = 0. Note that hab now is 
entirely two dimensional. Keeping this transverse metric in mind we can proceed in the 
same way as for the timelike case by constructing Bab = V^/ca. We quote below the equation 
for the expansion: 

^ + \^' + ^'-^' = -Rabk'^e (17) 

dX 2 

where the hatted quantitites are the expansion, rotation and shear for the null geodesic 
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congruence. The focusing theorem for null geodesic congruences follows in the same way 
as for timelike congruences, with the null convergence condition Rabk°'k^ > being the 
requirement. Using Einstein equations one can obtain the so-called Null Energy Condition 
Tabk^k'' > 0. Similar to the case for timelike congruences, we have corresponding equations 
for the evolution of shear and rotation for null geodesic congruences too. These are available 




E. The acceleration term and non afRne parametrisations 

The discussions above were exclusively for geodesic congruences. In the non-geodesic 
(i.e. timelike or null congruences) case, it is obvious that there will be differences. We 
state below, how the equation for the expansion changes via the addition of the so-called 
acceleleration term. The equation is now given as: 

^ + + Or^ - ^2 _ {^v'VbV'') = -RabV^'v' (18) 

where the fifth term on the R. H. S. is the acceleration term. Notice that this term is zero 
for geodesic congruences (zero acceleration). The average distance between the world lines 
in the flow is changed due to the divergence of the acceleration. The non-geodesic character 
of the flow also affects the equation for the shear. 

For a non-affine parametrisation of null geodesic congruences, with k^-Vak^ = f^k^ {k, 
a constant, defined through the above equation) the definition of the expansion changes : 
O = Vak"' — K and the Raychaudhuri equation for the expansion takes the form: 

^^^Q + l^Q^ + a'-u' = -Rabk^k' (19) 

The new feature here is the presence of the linear (in B) term. The conclusions on geodesic 
focusing however do not change, except for differences in the values of the expansion. 



F. A theorem for non rotating singularity free universes : Raychaudhuri's later 
papers 

Over the entire period of more than forty years, Raychaudhuri did not work much on 
the equations which bear his name today. Interestingly, he came back to have a look at 
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them once again in the late nineties with a series of papers [4l|, |42|. In these articles, he 
constructed a theorem on non-rotating singularity-free universes. The main inspiration 
behind these papers were the singularity-free cosmological solutions due to Senovilla j43| . 
which created a lot of interest and curiosity among relativists in the 1990's. In fact, the 
work of Raychaudhuri sets out to show that these solutions may not quite be physically 
relevant, though surely, mathematically correct. Let us now briefly recall the theorem of 
Raychaudhuri. 



The basic premise of j4l|, |42] is the use of spacetime averages of quantitites, defined as: 

"Ho i-k xMd'x 



ix) 



(20) 



rXQ rXi rX2 fX^^ \ I 74 

j-XQ j-xi J-X2 J-X3 y \9\(^ ^ 
Raychaudhuri shows that for a singularity-free non-rotating universe, open in all directions, 
the spacetime average of all stress-energy invariants, including the energy density vanishes. 
The proof is worked out using the spacetime averages of the scalars that appear in the 
Raychaudhuri equation for the expansion. Following the statement of the theorem, Ray- 
chaudhuri claimed that an observationally consistent universe cannot have a zero average 
density and hence one must necessarily give up the hope of having singularity-free solu- 
tions. Subsequent to Raychaudhuri 's work, Saa and Senovilla wrote a couple of comments 
which were mainly concerned with the converse of Raychaudhuri 's theorem and the ques- 
tion whether spacetime averages could really be a property which can distinguish between 
singular and non-singular models. In fact, one must note that the theorem does not say 
that if the spacetime averages are zero the spacetime must be non-singular. Therefore, 
the theorem cannot really be used to distinguish between singular and non-singular cos- 
mological models. Furthermore, Senovilla, in his comment, made a conjecture that spatial 
and not spacetime averages could be a distinguishing property between singular and non- 
singular models. More precisely, his claim was that for non-singular, non-rotating, globally 
hyperbolic and everywhere expanding models where SEC holds, the spatial averages of all 
stress-energy invariants must vanish. Therefore, if such spatial averages do not vanish then 
the model must be singular. More details and a recent proof of the conjecture is available in 



the article by Senovilla in this volume [4^ . Despite the conflict between whether spatial or 
spacetime averages was the crucial distinguishing factor between singular and non-singular 
models, it goes without saying that Raychaudhuri 's idea about averages being a deciding 
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factor will surely be remembered as a lasting contribution apart from his coveted equations. 



III. THE EQUATIONS IN DIFFERENT GEOMETRIES AND DIFFERENT THE- 
ORIES OF GRAVITY 

The Raychaudhuri equations discussed in the previous sections do not change as long as 
we respect the Riemannian (pseudo-Riemannian) metric structure of space or spacetime. 
However, in an alternative theory of gravity, the Einstein field equations are of course dif- 
ferent and hence the relation between the energy-momentum tensor and other geometric 
quantities do change. This results in a modification of the consequences that arise while 
analysing these equations. On the other hand, if the usual Riemannian structure of space- 
time changes, such as in the case of spaces with torsion (Einstein-Cartan-Sciama-Kibble 
theory), then the Raychaudhuri equations are surely different. We shall give an example 
of both these scenarios in the following two subsections. The equations in spacetimes with 
torsion is discussed in relatively greater detail because, this is one scenario where we actually 
notice a generalisation of the usual Raychaudhuri equations. 



A. Metric theories with symmetric connections 



As mentioned above, for theories with a symmetric connection, the Raychaudhuri equa- 
tions do not change, though the R.H.S. of the expansion equation when written using mat- 
ter stress-energy can be very different from what it is in GR. This change can surely affect 
geodesic focusing. It must be noted here that any change in the conclusions about geodesic 
focusing in this case is inherently due to the new solutions (spacetime geometries) of the 
modified Einstein equations for a given stress-energy. The abovementioned asp ects have 



47j as 



been analysed in the context of the Brans-Dicke and Hoyle-Narlikar theories 
well as the R + f3R^ theory 48|, on which we focus our attention below. 

The study of the Raychaudhuri equations in the presence of higher curvature terms in the 
action has been a subject of interest for a long time 49|]. The equation for the expansion, in a 
spacetime background solution of R + PR"^ gravity in the presence of matter (p oc a(t)~") has 
been analysed in ^|. The Strong Energy Condition (SEC) is examined and it is shown that 
the condition for a big-bang singularity changes in such a scenario. Recall that in GR, the 
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singularity theorem indicates the inevitable presence of singularities when the SEC (as well 
as some other conditions) are satisfied. However, at the microscopic scale, in the evolution of 
the expansion of a geodesic congruence, quantum effects are expected to show up - which, 
in turn, may completely change classical predictions. One such quantum effect leads to 
the appearance of higher curvature terms, namely, quadratic gravity 48|. Even though it 
might seem that the analysis in [so!] for quadratic gravity is classical, it is, in a broader 
sense, a semi-classical analysis. Classical solutions of quadratic gravity have been studied 
to explore the nature of sing ularities and it has been shown that the big-bang singularity 
may be avoided [511 • pOtl; quadratic gravity, i.e. R + [3R^ gravity, was considered as a 
a backreaction effect on pure Einstein's gravity . It was found that the SEC for R + (5R^ 
gravity is different from that of Einstein gravity. We briefly summarize the results of [sol 
below. 

In GR, the SEC follows from the use of Einstein's equations through the relation: 



RabV^v^ = 8nG 



T, 



ah 



1 



v^v" > 



(21) 



for all timelike v"". This, as mentioned before, leads to focussing. 

The author in [s^ is primarily interested in the cosmological singularity. To analyse the 
effects of the quadratic terms one only needs to flnd the new expression for the Ricci tensor 
Rab in terms of the energy momentum tensor Tab- Therefore, we require the modifled fleld 
equations for quadratic gravity. It is known that this equation is 



Gab + IGnGp ( ^R^gab - 2RRab 



'^R'ndab + '2R;a;b = SnGTab 



(22) 



In a perturbative analysis, we are primarily interested in the flrst order contribution from 
pR^ to the Raychaudhuri equation. After some algebra one arrives at. 



2f3G 



RabV^'v' = 
T^gabG — 2TTabG — T:n9ab — 2T-a] 



^Tgab) 



G{Tab 



(23) 



where G = StcG. 

Comparing the above equation with that in Einstein gravity, we flnd that an effective energy 
momentum tensor appears in the RHS of the above equation. Assuming a matter dominated 
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universe with a characteristic dependence on the scale factor (i.e., p = ^) and a local 
conservation of Tab leads to p = ^^f^- Thus, finally, we have. 



RabV^V 



n-2 ~ 



Gpn + I3G^ 3n(n - l)(n - A)Gpi + Qn\n - 4)A;a->„ 



+ 



{A' - 1) 



n ~ 



Gpn + PG' 2n\n - A)Gpi + 4n(n + 2)(n - 4)A;a~> 



(24) 



The above expression is the modified term which appears in the R. H. S. of the Raychaudhuri 
equation for the expansion. It is thereafter analysed in the context of different values of n and 
the possiblities of avoiding the big-bang singularity or its inevitable occurence are pointed 
out for different cases. 

We repeat once again that for any modified theory of gravity (eg. Brans-Dicke, Einstein- 
Gauss-Bonnet in higher dimensions, induced gravity, low energy effective stringy gravity etc. 
) as long the Riemannian structure of spacetime is respected, changes in conclusions related 
to geodesic focusing can arise only through the modified field equations and its use in the 
convergence condition. 

B. The Raychaudhuri equation in spacetimes with torsion 



The symmetric nature of the affine connection is one of the underlying assumptions of 
Riemannian geometry. This fact is also assumed while constructing GR. An antisymmetric 
connection may originate from the presence of spin matter fields in spacetime leading to a 



transition from the V4 to U4 manifolds 53] • This asymmetric part is known as torsion. 
Generalization of the theory of gravity in such a spacetime with torsion was proposed by 
Einstein-Cartan-Sciama-Kibble (ECSK).j53]. Absence of any experimental signature of 
torsion however is the primary criticism of such models although there has been a renewed 
theoretical interest in the context of superstring theories 5J] where spacetime torsion ap- 
pears in the form of a massless string mode. The third rank field strength of the massless 
second rank antisymmetric tensor field of string theory (known as the Kalb-Ramond field) is 
identified with spacetime torsion in the low energy limit of the theory. Torsion also appears 
naturally in a theory of gravity where twistors are used 5^ as well as in the supergravity 
scenario where torsion, curvature and matter fields are treated in an analogous way |5g] . It 



has been shown in several articles 



57l | that in string inspired models such a background with 
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torsion results in a departure from the experimentally predicted values of the well known 
phenomena like gravitational lensing, perihelion precession of planetary orbits, gravitational 
redshift, rotation of the plane of polarization of the distant galactic radiowaves etc. All 
the above arguments, and several more, compel us to include torsion in any comprehensive 
theory of gravity. In the context of geodesic congruences it has further been shown that the 
kinematical quantities- shear, rotation, acceleration, expansion and their evolution equa- 
tions are modified by the presence of torsion. This naturally leads to a generalization of 
the Raychaudhuri equations in presence of torsion leading to a more general understanding 
of the phenomenon of geodesic focusing. Here we give a brief review of the work presented 



m 



58l | where the role of torsion in modifying the Raychaudhuri equations and some of its 
implications have been discussed. 

The torsion tensor T^^'^ which is the antisymmetric part of the affine connection F^^, is 



'^ab" — n (^ab — Tba) = ^[ab] 5 (25) 



given by, 

1 
2 

where a, 6, c = 0, ... 3. 

In GR, T^f,'^ is postulated to be zero. 

From the torsion tensor one constructs the contortion tensor as, 

^ab'^ — ~'^ab'^ ~ '^'^ab + a = ~^a\ ? (26) 

This leads to a general expression for the affine connection given as. 



Kb = {lb} - K^b" , (27) 

where {^^} is the symmetric part of the connection (Christoffel symbols). With this modified 
connection the commutator of the covariant derivatives of a scalar field is 

V[,Vb]0 = -T„,'=V,0; (28) 

which is zero in absence of torsion. 

Similarly for a vector v'^ the commutator of the derivative gives, 

V[aVb]v' = Rabdv" - 2T,/V,t;^ (29) 
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where the Riemann tensor is defined as, 



p d p) prf f) 'pd I 'pd -pe -pd -pe 

^abc ~ ^a-*^ be '^b'- ac ae-*- be be'- ae- 



(30) 



The contribution of torsion, to the Riemann tensor, is exphcitly given through the following 
expression: 

Rabe' = Rabe\{}) " ^ aK,, ' + VbKj + KJK,; - KJKJ (31) 

where Rab,f'{{}) is the tensor generated by the Christoffel symbols. The symbols V and V 
have been used to indicate the covariant derivative with and without torsion respectively. 
From Eq. (!3T]) . the expressions for the Ricci tensor and the curvature scalar are 



Rab = RabiU) - 2V„T, + V,KJ + KJK,^^ - 2nK, 



(32) 



and 



where, 



R = RiU) - 4V,T'^ + K,,,K'''' - ATaT\ 



Ta can be either time-like, space-like or light-like. 



1. Contributions of torsion to shear, expansion, vorticity and acceleration 



(33) 



(34) 



6l| , 62| . In 631] an infiationary 



6J] and [65| addressed 



Beginning with the behaviour of fiuids and moving on to the initial singularity problem in 
cosmo logical models, the Raychaudhuri equations have been shown to be the key equation to 
explore the role of torsion in such diverse phenomena, jsol . 60], 
model of the universe in the context of ECSK theory was considered, 
the crucial role of Raychaudhuri equation in the context of a gauge invariant formalism 
for cosmological perturbations in theories with torsion. We now explain the kinematical 
quantities mentioned before, in a scenario with torsion. 

Torsion in a space-time modifies the definition of kinematical quantities. The covariant 
derivative of the four velocity Ua can be decomposed as. 



&ab + -^habQ + UJab " f^a«fe 



(35) 
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where hat = Qab + UaUb and 

e = Vat/° = e - 2T'^f/e, (36) 

= KhiVi^JJd) = + 1KKK^^^{V,, (37) 

(ba^> = KhtVlJJd] = tOab + 2hlhtK^^^{Ue, (38) 

and the acceleration 

ae = U'^WaUc = a, + U^KjUd. (39) 

The quantities without the tilde are the value of the corresponding expressions in a spacetime 
without torsion. 

2. The Raychaudhuri equation 

Using the identity for the four- velocity Ua {UaU"' = — 1), 

U'V.VbUa = V,{U'VbUa) - VM'VbUa (40) 

and from Eq. fl30l) 

U'V.VtUa = U'VbVJJa + R^JUdU' - 2UX, ^Vrft/, (41) 
we find the equation 

(l ~ 2~ ^ 2~ ^ 

- I g ^ca© + -00-ca + -QuJca + 2a^UJba 

- R^baUdU^ - 2U%b'' [IhdcQ + + Codc - Uda^j . (42) 

Contracting the indices in Eq. fl42l) . one obtains a general expression for the Raychaudhuri 
equation for the expansion, in the presence of torsion as, 



+ a^, + - Uda'^) (43) 
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This is the most general form of Raychaudhuri equation in presence of torsion. Simpler 
versions of this equation have been discussed in 61], 66], [s^jjo^]. It is obvious that there 



would be corresponding equations for shear and rotation too, which we do not mention here. 
Interestingly, it can be shown that if we have torsion as a phantom field (a scalar field with 
a negative kinetic energy term) through the dual of the third rank tensor field strength of a 
string inspired Kalb-Ramond field, one can get a positive contribution in the R.H.S. of the 
above equation, which, in turn, may be useful in eliminating singularities. 

IV. THE EQUATIONS IN DIVERSE CONTEXTS 
A. General Relativity and Relativistic Astrophysics 

We have already discussed one of the main uses of the Raychaudhuri equations-the 
geodesic focusing theorem. Though it is an entirely geometric result, its use is mostly 
confined to the domain of GR. It is also quite obvious that the null version of these equa- 
tions are useful in the study of gravitational lensing. The focal point in this case is nothing 
but the intersection of trajectories representing light rays and is known as the caustic of 
the bundle of trajectories. The optical scalars (Sachs scalars) are the quantitites of interest 
and their evaluation enables us to understand the nature of the null geodesic flow (light ray 
bundles). Is there anything else one can say apart from focusing, which will be relevant 
for lensing. Some of these issues have been addressed in 68|. For example, it is possible 
to rewrite the Raychaudhuri equation for null geodesies in a form involving the angular 
diameter distance cIa- 

% - - V («' 

where ^Inrf^ = |B. Thus, for the case of zero shear one can determine dA entirely in 
terms of geometry. In [GSj the authors have asked the question 'how do voids affect light 
propagation?' and made use of the above relation to provide an answer. 

A second example within relativistic astrophysics, is a study of crack formation 0, [t^ 
using this equation. This is an interesting application, quite unique- but not pursued much 
later. The basic question here is - when will a spherical object develop cracks? Appearance 
of total radial forces of different signs in different regions in a perturbed configuration leads 
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to the occurrence of such cracking, which, in turn leads to local anisotropy of fluid/emission 
of incoherent radiation. Using the Einstein and Raychaudhuri equations it is possible to 
express the net radial force R (after perturbation) in terms of pj, p, gij. In this way a 
criterion for cracking can be obtained. 



In 1995, Jacobson 7l|] in a rather unique paper demonstrated how one might view the 
Einstein field equation as a thermodynamic equation of state. In this calculation, Jacobson 
starts out by arguing that energy flux across a causal horizon is some kind of heat flow 
and the entropy of the system beyond is proportional to the area of the horizon. The heat 
flux 6Q is related to the stress-energy tensor whereas the area variation is related to the 
expansion of a bundle of null geodesies. He then makes use of the Raychaudhuri equation 
in its linearised form (ignoring the term) to write down the expansion as an integral 
over the Rijk^y . Thereafter, with the input of the entropy-area relation he arrives at the 
Einstein equation as an equation of state. It is important to note here that Jacobson uses 
the geometric content of the Raychaudhuri equation and views it as more fundamental than 
the Einstein equation, which, in his approach is a derived relation. Extensions of Jacobson's 
work in the setting of non-equilibrium thermodynamics have appeared recently [72]. 

Scattered around the literature, are innumerable instances where the equations have 
been used in the context of GR and astrophysics. Prominent among them is its use in 
black hole physics-in studying the properties of black holes and in deriving the laws of 



black hole mechanics (see ISj and references therein). Apart from black holes, we mention, 
en passant, a few other situations which we thought were interesting: (i) use in a fiuid- 
fiow description of density irregularities in cosmology [73j], (ii) quantum gravitational optics 
and an effective Raychaudhuri equation [t^, (iii) magnetic tension and gravitational collapse 
[75I , (iv) the effective Einstein and Raychaudhuri equations derived from higher dimensional 
warped braneworld models [tG]. 

B. The Capovilla Guven equations for relativistic membranes 



An obvious question, which was never asked till the work of Capovilla and Guven [77 1 
appeared in the scene is: what happens if we consider a congruence of extremal surfaces 
instead of geodesies (extremal curves)? Are there similar Raychaudhuri equations? 

To address this issue, we must first find out how to generalise the notions of expan- 
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sion, rotation and shear for the case of a family of surfaces. Unhke a curve, a surface is 
parametrised by more than one parameter. Thus, it is natural to imagine an expansion, a 
rotation and a shear along each of these independent directions. Introducing a separate label 
for the surface coordinates, wc find that we now have G", S"^- and Q^j (here denotes 
the indices representing the normals to the surface). The relation between the i,j indices 
and the spacetime indices (say /i, u) is given through the embedding of the surface. Let us 
now make things more concrete and explicit. 

Define a D-dimensional surface in a N dimensional background through an embedding 
J.H _ x^{^"'). Ej^ constitute the tangent vector basis chosen such that g^^E^E^^ — r]ab (a,b 
run from 1 to D). n^"^ arc the normals, with g^yU^'^n'^^ = 5''^ (i,j run from 1 to N — D). Also 
g^yii^^E'^'^ = 0. K"-^^ are the N — D extrinsic curvatures (one along each normal direction) 
The embedded surface is minimal provided Tr{Ki) — 0. 

With the above definitions, one can follow the derivation for the case of geodesic curves 
and obtain the equation for the generalised expansion (assuming zero values for the Ejj and 

v«e„ + ^^^e„e" + (m^)| = o (45) 

where 

(m^) i = K'^'^Ka^ + R^^.aE'^^EX'nt (46) 

Note that the above equation is a partial differential equation-the reason being that we 
require more than one parameter to describe a surface. Structurally, the equation is similar to 
the original Raychaudhuri equation for the expansion - one can easily note this by comparing 

the two. However, there are crucial differences - one of which involves the appearance of 
extrinsic curvature terms. It is possible to rewrite this equation in a second order form by 
choosing = Va-F. We obtain: 

VaV"F+ (m2)^F = (47) 

which, resembles a variable-mass wave equation on a D-dimensional surface with a non- 
trivial induced line element. For special cases, one can work out focusing criteria, although 
the notion of focusing will be largely different for surface congruences. For the case of string 
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world-sheets, it is possible to make use of the conformal character of any two-dimensional 
line element and rewrite the equation as: 



g-|^ + tf(r.-)(M=)> = (48) 
where r, a are the coordinates on the string world sheet and the induced metric is ds^ = 



It can be shown 



78| that the criterion for focusing is given as (where is the Ricci 



scalar of the worldsheet) 



-^R + R^^E^'^E"^ > (49) 



which is a requirement for the function F to have zeros. However, there are many questions 
related to focusing of surfaces which remain un-answered. Some of these have been recently 
addressed in [tQ^ . Earlier references on examples of solutions of the Raychaudhuri equations 
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80| 



for surface congruences and related issues are available in 

It must be mentioned here that the above generalised equation for the expansion vector 
field is a subset of a full set of equations which include those for the generalised shear and 
rotation too. Moreover, the analysis is entirely for extremal membranes of Nambu-Goto 
type. Here too, if the action changes, we find that the equations change-an example is 
available in j81|. 

C. Quantum field theory 

In recent times, the kinematic quantitites (expansion, shear and rotation), as well as the 
Raychaudhuri equations, have appeared, quite unexpectedly, in the context of quantum field 
theory. We outline below, some of these scenarios briefly. 

1. The Langevin-Raychaudhuri equation 

A couple of years ago, Borgman and Ford investigated gravitational effects of quantum 



stress tensor fluctuations |82|, [83(1 • They showed that these fluctuations produce fluctua- 
tions in the focusing of a bundle of geodesies. An explicit calculation using the Raychaudhuri 
equation, treated as a Langevin equation (ignoring the 0^ term by a smallness assumption 
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converts the Raychaudhuri equation to the form = /(A), a Langevin-type equation) was 
performed to estimate angular blurring and luminosity fluctuations of the images of distant 
sources. The stress-tensor fluctutations were obtained assuming the case of a massless, 
minimally coupled scalar field in a fiat background in a thermal state. Scalar field fiuctu- 
ations drove the Ricci tensor fiuctuations (via the semiclassical Einstein equations), which, 
in turn led to fiuctuations in the expansion 0. These authors also made some numerical 
estimates for the quantity ^ (the fractional luminosity fiuctuation) and pointed out pos- 
sible astrophysical situations (gamma-ray bursts, for example), where an observable value 
of this quantity might exist. However, it is true that their work has many assumptions 
(ignoring shear, rotation as well as the 0^ term etc.) and much further analysis is required 
to have a clearer picture of the possibility iof actually seeing such effects in the real universe. 
Nevertheless, the idea, on the whole, is novel and interesting in its own right and provides 
us with another application of the Raychaudhuri equation in a very different scenario. 

2. RG flows in theory (coupling) space 

Till now we have been looking at geodesic flows in spacetime. However, one might also 
contemplate geodesic flows in flctitious spaces where a metric is deflned. Such spaces may 
correspond to certain physical scenarios. We highlight one example here in order to show 
how ubiquitous the notions of expansion, rotation, shear and the Raychaudhuri equation 
are. 

For a given Lagrangian fleld theory, the set of couplings can form a space which is known 
as theory/coupling space. A curve in such a space will therefore represent a flow of couplings. 
In such a space, we can deflne a distance function using two-point functions integrated over 
physical spacetime. Such a deflnition of metric goes back to the well-known Cramer-Rao 
metric in probability theory and has also been highlighted in the context of fleld theory by 
Zamolodchikov as well as O'Connor-Stephens (see [8^ and references therein). 

Once we have a metric, we can use it to deflne derivatives of a vector fleld. The vector 
fleld of interest here is the so-called (3 — function vector fleld, which generates the RG flow. 
The covariant derivative of this vector fleld when split into the usual trace, anti-symmetric 
and symmetric traceless parts deflne the usual expansion, rotation and shear. However, we 
also know that the /^-function vector fleld must satisfy the conformal Killing condition (this 
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is the geoemtric form of the RG equation). These facts together imply the result that the 
focal point of the congruence generated by the /3-function vector field must necessarily be a 
fixed point (i.e. all Pa ('a' labels the coupling space coordinates) vanish there) [84 1 

This line of thought is useful in obtaining generic results in field theory -the under- 
standing of shear and rotation in the context of RG flows might provide a better geometric 
understanding of RG flows. 



3. Holography, c-function and the Raychaudhuri equation for the expansion 

The holographic principle has recently played a crucial role in our understanding of 
quantum aspects of gravity. The principle states that the information of gravity degrees 
of freedom in a D dimensional volume is encoded in a quantum field theory defined on the 
(D — l) dimensional boundary of this volume. As an immediate realization of this, it is shown 
that the Renormalization Group flow equation for the /3-function of a four dimensional 
quantum gauge field theory defined on the boundary of a five dimensional volume can 
be described by geodesic congruences in a scalar-coupled five dimensional gravitational 
theory. Such a gauge-gravity duality was proposed in a more general framework through the 
Maldacena conjecture and can be elegantly described through the Renormalization Group 
equation with the bulk coordinate (or the holographic coordinate) as the Renormalization 
Group parameter. It is shown that if the central charge or the c-function in a quantum 
field theory evolves monotonically under RG then the holographic principle indicates that 
in the corresponding dual gravity in five dimensions, the picture is realized through the 



Raychaudhuri equation governing the monotonic f 



ow o 



the expansion parameter B for 



85|, |86||87|. The central charge of the 



the geodesic congruences in the gravity sector 
boundary theory, or the c-function, is a measure of the degrees of freedom of the theory. 
As a consequence it is also a measure of the entropy of the black hole in the dual gravity 
sector which is related to the area of the horizon through the Bekenstein-Hawking formula. 
The effective central charge is a function of the couplings of the theory, which monotonically 
decreases as one fiows to lower energies through the RG equation. The fixed point described 
by the boundary conformal field theory corresponds to the extrema of this function. It turns 
out that the null geodesic congruences can be used as a probe to decode the holographically 
encoded informations. 
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Consider a D-dimensional spacetime with a negative cosmological constant where the 
metric gab is fohated by appropriate choices of constant time surfaces. If we now focus on a 
(D — 2) dimensional spacetime surface M at a fixed time t, then one may construct a null 
vector field on the light sheets (consisting of spacetime points scanned by null geodesic 
congruences) which is orthogonal to M such that m"'na = — 1, where Ua are tangents to the 
geodesic. The metric on M is given by 

hab = 9ab + n^arrib) (50) 

Defining, 

Bab = DbUa (51) 
such that Bab = B^a (Frobenius theorem ensures the symmetry property) and 

e = TrB (52) 

it can be shown from the Raychaudhuri equation for the expansion that the geodesies con- 
verge since, 

e < (53) 

The question that arises now is - how is the information on the light sheet encoded in M? 
Assuming that the spacetime admits a timelike Killing vector field, we may visualize the time 
flow of M to constitute the {D — 2) + 1 — [D — 1) dimensional boundary of the D dimensional 
bulk. The c-function of the corresponding dual gauge theory sits on the boundary. The RG 
flow to lower energy scales is then associated to the motion along the converging congruences 
of null geodesies described by Raychaudhuri equation. In this sense, the null geodesies act 
as a probe for the dual theory at a lower energy scale. Pictorially, it may be described 
as follows. The lower the energy scale, the RG takes us into the gauge theory, deeper 
inside the bulk we move along null geodesic congruences following Raychaudhuri equation. 
The caustic point where the geodesies meet correspond to the fixed point of the RG flow. 
This establishes a remarkable significance of the Raychaudhuri equation in describing the 
holographic principle. 

As described previously the c-function tells us about the degree of freedom of the system, 
which, in turn, is related to black hole entropy through the area law. Therefore, the c- 
function is naturally related to the area function A{r) as, 

c(r)=^(r)/4 (54) 
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In the context of the recently developed attractor mechanism 88|]in describing black holes 
in a scalar coupled gravitational theory, one uses this flow of the c-function to decode some 
interesting properties of both supersymmetric and non-supersymmetric attractors. The 
Raychaudhuri equation once again plays a crucial role. It is shown that in any spherically 
symmetric scalar coupled static asymptotically flat solution, c(r) decreases monotonically 
as one moves radially inward from inflnity. It is further shown that the minimum value of 
c(r) corresponds to the entropy of the horizon 89l| . In order to establish the c-theorem in 
this context once again the Raychaudhuri equation is used. Taking a congruence of null 
geodesies, the expansion is given by, 

e = ^ (55) 

(where A is the affine parameter). Using the null energy condition, 

Tabk^k^ > (56) 

and the Raychaudhuri equation, we get 

The c-theorem can thus be proved on very general grounds. 

In summary, these rather novel applications of the equations in the context of quantum 
fleld theory seems to assert once again the generality and the ubiquitous nature of the 
equations. 



V. SUMMARY AND SCOPE 



We have reviewed, in this article, the Raychaudhuri equations and its applications in 
diverse contexts. It must be admitted that since its inception, the equations have used 
extensively to enhance our understanding of various situations within, as well as outside the 
scope of GR. The primary reason behind its wide-ranging applicability is (as emphasised 
several times in this article) the fact that the equations encode geometric statements about 
flows. Since flows appear in many different contexts in physics, it goes without saying, that 
the equations will be useful in furthering our understanding in different ways. 

Among the applications within classical GR, the utility of the equations in understanding 
the singularity problem, is surely the most prominent one. In astrophysics, we have quoted 

29 



the use of the equations in the context of lensing, cracking of self-gravitating objects. At a 
more fundamental level, Jacobson's work exploits the equation and, in some sense, makes 
use of its geometric nature to arrive at a new way of looking at the Einstein equations. 

On the other hand, in a non-Riemannian spacetime, we have seen how the equations 
change -in particular, in spacetimes with torsion. Torsion appears in disguise in string 
theory and is therefore not totally irrelevant today, even though observations might have 
ruled out the Einstein-Cartan-Sciama-Kibble theory long ago. 

Furthermore, what has come as a pleasant surprise in recent times, is the utility of these 
equations (and the quantities that appear in the equations: expansion, shear and rotation) 
in the context of quantum field theory. In addition, the generalisation of these equations to 
surface congruences is also an interesting development, despite the fact that not much has 
been done on such generalised Raychaudhuri (Capovilla-Guven) equations. 

Despite the large body of work on or based on the Raychaudhuri equations, there still 
remain many un-answered questions. In conclusion, we list a few of them. The choice is 
surely biased by our own perspective. 

(a) Normally, in fluid mechanics, the independent variable is the velocity field. The 
Raychaudhuri equations arc for the gradient of the velocity ficld-thcy are one derivative 
higher. We mentioned before that they are essentially identities and become equations 
when we choose a geometric property (eg. vacuum, Einstein space etc.) or use Einstein's 
field equations of GR. Is it possible to set-up and solve the initial value problem for this 
coupled system of equations such that we know the characteristics of a geodesic flow in a 
given spacetime geometry once we specify the initial conditions on expansion, shear and 
rotation? Similar analysis is done while deriving the focusing theorem but can we do it by 
analysing the full system of equations analytically/numerically. Following this approach one 
might be able to reconstruct the gradient of the velocity field and also obtain the velocity 
field itself. 

(b) In the context of RG fiows, is it possible to understand the role of shear and rotation 
and use it to improve our geometric analysis of such fiows? 

(c) For surface congruences, can we obtain equations for families of null surfaces? Or, for 
surfaces defined by the extremality of actions other than Nambu-Goto (rigidity corrections, 
Willmore functionals etc.) 

(d) In quantum mechanics, we talk about the fiow of probability. Imagining this as a fiuid. 
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can we obtain expansion, rotation and shear of such flows and then write down conclusions 
based on the corresponding Raychaudhuri equations? 

(e) Electrodynamics, is, after all, the physics of the electric and magnetic vector fields. 
Suppose we construct the gradient of the electric/magnetic fields, split it up into expansion, 
rotation and shear and write down Raychaudhuri equations for the so-called Faraday lines 
of force! What does such an analysis tell us? 

To sum up, we would like to conjecture that wherever there are vector fields describing 
a physical/geometrical quantity, there must be corresponding Raychaudhuri equations. We 
believe, there are many avatars of the same equations, of which, till today, we are aware of 
only a few. 
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